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Abstract 



The properties of three-jet events with total transverse energy greater than 320 GeV 
and individual jet energy greater than 20 GeV have been analyzed and compared 
to absolute predictions from a next-to-leading order (NLO) perturbative QCD cal- 
culation. These data, of integrated luminosity 86 pb -1 , were recorded by the CDF 
Experiment for pp collisions at ^/s = 1.8 TeV. This study tests a model of higher 
order QCD processes that result in gluon emission and can be used to estimate the 
magnitude of the contribution of processes higher than NLO. The total cross sec- 
tion is measured to be 466 ± 3(stat.)lyo ? (syst.) pb. The differential cross section is 
furthermore measured for all kinematically accessible regions of the Dalitz plane, in- 
cluding those for which the theoretical prediction is unreliable. While the measured 
cross section is consistent with the theoretical prediction in magnitude, the two differ 
somewhat in shape in the Dalitz plane. 
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In perturbative QCD, hard scattering of the constituent partons in the proton 
and antiproton results in events with large total transverse energy, J^Et- Outgoing 
scattered partons hadronize and may be detected as hadronic jets. Three-jet events 
can be produced when a hard gluon is radiated from any of the initial, intermediate, 
or final state partons in an event with two primary outgoing partons. 

We analyze here some properties of the cross section for three-jet event production 
in proton-antiproton collisions at the Fermilab Tevatron Collider at center-of-mass en- 
ergy 1.8 TeV. The data, which were recorded by the Collider Detector at Fermilab 
(CDF) [1], are compared with predictions by the first complete next-to- leading or- 
der (NLO) QCD generator, Trirad [2], for hadronic three jet production at hadron 
colliders. We compare the measured and predicted absolute cross sections to test 
our understanding of the higher order QCD processes that result in gluon emission 
and to estimate the magnitude of the contribution of processes higher than NLO. 1 
In some kinematical regions, we provide a measurement of the cross section where 
the theoretical prediction is not reliable; this measurement may be a useful guide for 
theoretical calculations. We also compared the shapes of the measured and predicted 
cross sections when normalized, to examine the sensitivity of the cross section to 
variations in the value of the strong coupling, a s . The data sample corresponds to an 
integrated luminosity of 86 pb _1 collected during the 1994-1995 run (Run lb). 

A previous paper [4] examined a smaller dataset and was limited to a comparison 
with leading order theoretical calculations [5]. A subsequent analysis [6] compared 
a larger dataset to predictions from the HERWIG [7] parton shower Monte Carlo 
program and to the NJETS [8] leading order 2 — > N parton-level prediction. The 
NLO calculation used here has the benefit of reduced renormalization scale depen- 
dence (and consequently lower systematic uncertainty) as well as a more reliable 
description of multijet production throughout phase space. This study expands upon 
the previous investigations by comparing the data to absolute cross section predic- 
tions. The measurements presented here include differential cross sections that may 
be useful constraints upon parton distribution functions. 

We use a coordinate system with the z axis along the proton beam, transverse 
coordinate perpendicular to the beam, azimuthal angle 0, polar angle 9, and pseudo- 
rapidity i] = — lntan(0/2). The analysis uses the CDF calorimeters [9], which cover 
the pseudorapidity range \rj\ < 4.2. The calorimeters are constructed in a tower ge- 
ometry and are segmented in depth into electromagnetic and hadronic components. 
The calorimeter towers are 0.1 unit wide in rj. The tower widths in are 15° in the 
central region and 5° for \r]\ greater than approximately 1.2. 

We begin by considering events from the data sample selected by the trigger 
requirement > 175 GeV. We refer to this 175 GeV as below. Event recon- 

struction uses a cone algorithm [4] described in more detail below. The transverse 

1 While no quantitative estimate of the contribution of next-to-next-to-leading order processes to 
the cross section is available at this time, considerable progress has recently been made in calculating 
two loop 2^2 parton processes [3], important groundwork for the future. 
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energy is defined as £r = E sin 6, where E is the scalar sum of energy deposited in 
the calorimeter within a particular cone and 8 is the angle between the beam direction 
in the laboratory frame and the cone axis. All calorimeter energy clusters [4] in the 
event with > 10 GeV are summed. The three leading jets in the laboratory frame 
are used as the basis of transformation into the three-jet rest frame. In the three-jet 
rest frame, the incoming partons are, by convention [10], labeled partons 1 and 2, and 
their momenta are designated pi and p 2 , respectively. The highest energy jets in this 
frame have energies labeled E 3 , £4, and E 5 and are ordered such that E 3 > E A > E 5 . 
The outgoing partons associated with these jets are correspondingly labeled partons 
3, 4, and 5. 

A three-jet system in the massless parton approximation can be uniquely described 
by five independent variables (see Figure 5 in [11]). We use the following: 

1. the invariant mass of the three-jet system, m 3 j 

2. the cosine of the angle #3 between the average beam direction (pay = Pi — P2) 
and parton 3 in the three-jet rest frame: 

a* Pav • P3 
cos 6> 3 = 



IpavI IpsI 

3. the cosine of the angle ip* between the plane containing the average beam direc- 
tion and parton 3 and the plane containing partons 3, 4, and 5 in their center 
of mass frame: 

, * (P3 x Pay) ■ (P4 x Ps) 

COS 1p = ——; — ^— 

\P3 X PAY] \P4 X p 5 \ 

4. the Dalitz variable X 3 (see below) for the leading jet, and 

5. the Dalitz variable X 4 (see below) for the next-to-leading jet. 

The invariant m 3J is calculated by sorting jets by their energies in the laboratory 
frame, boosting to the rest frame of those with the three highest energies, re-sorting 
jets by energy in that frame, then computing m 3 j = X^=3 E{, where the E\ are the 
energies of jets 3, 4, and 5 in the rest frame. We have investigated the probability 
that a jet with energy less than the weakest of the three jets in the laboratory frame 
may have an energy greater than E 5 in the 3-jet rest frame from which it is excluded 
by this algorithm. The restriction imposed by the cut on full trigger efficiency (see 
below) makes this probability negligible. 

The Dalitz variables, X- x , are defined as X\ = 2 • E[/m 3 j, (i = 3, 4, 5). Momentum 
conservation restricts the ranges of the Dalitz variables to 

I < A', < 1, 

- < X4 < 1, and 

< X 5 < -. 
- 5 - 3 
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A set of trigger and offline requirements [12] rejects events associated with cosmic 
rays, beam halo, and calorimeter malfunctions. Events are required to have a recon- 
structed primary vertex, defined as the vertex with the largest J2i P\ (where P ; is the 
total momentum of particle i leaving the vertex in the event), within \z\ < 60 cm. 
Events are defined to have resolved multiple interactions if a second vertex with at 
least ten associated tracks is reconstructed in the vertex track detector, and if that 
vertex is separated from the primary one by at least 10 cm. Because multiple interac- 
tions can change the jet multiplicity in an event, for example, misidentifying two-jet 
events as three-jet events, events with resolved multiple interactions are removed. 
The number of events with unresolved multiple interactions in which an additional 
jet could be misidentified is estimated to be less than 2% so no correction for them 
is applied. The resulting effective total integrated luminosity of the data sample is 
77±4 pb _1 , where the uncertainty reflects both the overall luminosity uncertainty 
(4.2%) and the uncertainty (0.5%) associated with the removal of resolved multiple 
interactions. 

An iterative cone algorithm [4] with cone radius R = \J (Ar/) 2 + (A</>) 2 = 0.7 is 
used to identify jets. Here Ar] = i] 2 — i]i and A<p = (j) 2 — <j>i- The subscripts 1 and 
2 correspond to the axes of the cone and calorimeter tower, respectively. Jets that 
share towers are combined if the total of the shared towers is greater than 75% 
of the Ei of either jet; otherwise the towers are assigned to the nearest jet. Jet 
energies are corrected [4] for errors in the absolute and relative energy scales and 
for additional energy associated with the underlying event. Since partons that are 
radiated out of the cone lead to the same losses in the theoretical calculation and 
in the data, out-of-cone corrections are not applied. The E T of a jet is calculated 
from the reconstructed position of the primary event vertex. All three leading jets 
are required to have E^ > 20 GeV and < 2.0. Events with fewer than three jets 
are rejected. To avoid collinear soft gluon instability in the iterative jet clustering 
algorithm [13], a cone overlap cut is imposed: events are rejected if the distance AR 
in rj-cf) space between the axes of any two of the three leading jets is less than 1.0 (see 
Figure 5 of [4], which shows that this selection requirement reduces to approximately 
zero the probability of the two jets being merged by the clustering algorithm). To 
exclude regions in which the geometrical acceptance [11] is less than about 95%, we 

require |cos 0|| < y^l — Je$£ / c 2 / m 2 j) , where m 2 j is the mass of the two leading jets 
in the three-jet system and is defined analogously to m 3 j. 

We require full trigger efficiency, which occurs when ^3 jets Et > 320 GeV, where 
the sum is over the three highest energy jets in the event with corrected E T > 
20 GeV [14]. The data are compared to the theoretical prediction by sorting events 
into bins of size 0.02 x 0.02 in X3-X4 space, the Dalitz plane. Figure 1 shows the 
Dalitz distribution of data that remain after all of the selection requirements have 
been applied. 

Before the final binning is done, the data are corrected for the effects of the 
combination of detector resolution and energy mismeasurement. A correction factor 
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Figure 1: The three-jet data, after all selection requirements have been applied. The 
energy correction procedure (see text) has not been applied. The figures at the corners 
of the distribution represent typical three-jet topologies in those regions of the Dalitz 
plane. 
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is determined for each bin in the plane as follows. A sample of events is generated 
at the parton level with the HERWIG Monte Carlo. The final state partons are 
hadronized. The events are then binned in the Dalitz plane. The same events are 
next passed through the CDF detector simulation and rebinned. For each bin the 
ratio of the number of events after and before detector simulation is computed. This 
ratio (ranging from 0.85 to 1.5) is the factor subsequently used to correct the number 
of events in each data bin. The data are also corrected for the z-vertex cut efficiency 
and then normalized to the effective total luminosity. 

The principal sources of systematic uncertainty [15] on the cross section are those 
on the absolute and relative (77- dependent) jet energy scales. The uncertainty on 
absolute jet energy derives from the resolution on the calibration of the calorimeter 
(uncertainty 1.3%-1.8%, and .Bp-dependent), the uncertainty associated with choice 
of jet fragmentation model (decreasing from 1.7% to 1.2% with increasing E T ), the 
uncertainty associated with calorimeter stability over time (1%), and the uncertainty 
on the correction for the contribution of the underlying event (1 GeV). The uncer- 
tainty on the relative jet energy scale ranges from 2% to 6%. Uncertainties are also 
associated with the measurement of the effective total integrated luminosity (4.2%) 
and with the z- vertex cut efficiency (2%). There is also an uncertainty of less than 
5% associated with the implementation of simulated events in the correction proce- 
dure. The upper (lower) limits on these uncertainties are added (subtracted) from the 
four-momenta of the jets in the data sample to obtain the systematic uncertainties 
on the cross section associated with each contribution. The uncertainties are then 
combined to produce the total experimental systematic uncertainty reported for each 
bin in Tables 1 to 14 below. 

The Trirad calculation consists of 2 — > 3 parton processes at one loop and 2^4 
parton processes at tree level. For gluons g, incoming quarks q, and outgoing quarks 
Q or Q', the subprocesses involved are gg — > ggg, qq — > ggg, qq — > QQg, and those 
related by crossing symmetry, all computed to one loop; and gg — > gggg, qq — > gggg, 
qq — > QQgg, and qq — > QQQ Q' and the crossed processes computed at tree level. 
The program uses the "subtraction improved" [13] phase space slicing method to 
implement infrared cancellation. 

The cross section is predicted with the CTEQ4M [16] parton distribution function 
(PDF) for each bin in the Dalitz plane. The result is multiplied by the effective 
total integrated luminosity of the data to predict a number of events in each bin. 
We restrict the prediction to bins for which X 3 < 0.98; this is necessary as the 
perturbative expansion is not reliable where the three-jet configuration approaches a 
two-jet configuration. The comparison between the data and the calculation is made 
for 215 bins. 

Figure 2 compares Dalitz distributions of the data and the absolute theoretical 
prediction. The theoretical distribution is more strongly peaked — a trend that persists 
in comparisons with all members of the CTEQ4 family 2 of PDFs. This trend, in which 

2 The CTEQ4 family includes CTEQ4A1, CTEQ4A2, CTEQ4M, CTEQ4A5, and CTEQ4A6, 
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CDF Run 1 Data 




Figure 2: The event density in the Dalitz plane for (a) the data and for (b) the predic- 
tion by the NLO Monte Carlo calculation with CTEQ4M, normalized to luminosity. 

the edges of the Dalitz plane are more populated by data than by the prediction, may 
give some indication of the size of the higher order contributions to the cross section. 

The data and theory are compared in two different ways. In Figure 3, we compare 
the shapes of their Dalitz distributions by normalizing the data and theory predictions 
to the same number of events. In Figure 4, we normalize theory to the experimental 
luminosity and compare the absolute values of the cross sections that are observed 
and predicted. In both figures, the prediction is made using the CTEQ4M parton 
distribution function, and the difference between observed and predicted number of 
events, scaled by the number of predicted events, is computed. 

The theoretical prediction for the cross section, using CTEQ4M and all bins in 
the Dalitz plane but those with X 3 > 0.98, is 473 ± 2(stat.)l!(scale)t|8( PDF ) P b - 
The theoretical uncertainty associated with choice of renormalization and factoriza- 
tion scales, /iR and /xf respectively, is estimated by varying the scales, whose default 
value is Et, to values of Et/2 and 2Et while keeping /xr, = /xp. The theoretical 
uncertainty associated with choice of PDF is estimated from the spread in the pre- 
dictions generated with all members of the CTEQ4 family. The measurement is not 
sensitive to the value of a s as is also shown in [17]. The measured cross section, using 
all bins in the Dalitz plane but those with X3 > 0.98, is 458±3(stat.)ljg 3 (syst.) pb. 
This is consistent with the theoretical prediction and with a previous CDF measure- 
ment [14] after corrections are made for the efficiencies of additional cuts introduced 

which differ in the value of a s input to their global fit, and CTEQ4HJ, for which a higher statistical 
emphasis was given to the high Et data from CDF. 
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Figure 3: The fractional difference between the data and the theoretical prediction, 
using CTEQ4M, as a function of X^, averaged over X3. The vertical bands show the 
systematic uncertainties. The error bars show the statistical uncertainties in cases 
where those are larger than the size of the symbol used. The prediction is normalized 
to the data to facilitate a comparison of the shapes of the distributions. 
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Figure 4: The fractional difference between corrected data and the NLO prediction, 
using the CTEQ4M parton distribution function, as a function of X3 for various X4 
bins. Error bars reflect statistical uncertainty for cases in which it is larger than 
the size of the symbol used. Shaded bands indicate systematic uncertainty. The 
prediction is normalized to the luminosity of the data. 
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in this analysis. The measured cross section, using all bins in the Dalitz plane, is 
466±3(stat.)t^ 7 (syst.) pb. 

Tables 1-14 summarize the measured cross section for every kinematically allowed 
bin in the Dalitz plane, including those for X 3 > 0.98. The measurements at high 
X 3 may provide useful constraints on future theoretical models in that region. The 
absolute predicted cross sections, using CTEQ4M and renormalization and factoriza- 
tion scales \i = E?, are also provided for bins with X 3 < 0.98. In a few bins, the 
predicted value is extremely small and dominated by theoretical uncertainties; for 
these bins no predicted value is quoted. The uncertainties on the measured values 
are the quadrature sum of statistical and systematic uncertainties; the uncertainties 
on the prediction reflect the statistics of the Monte Carlo sample. 



x 3 


x A 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.67 


0.67 


5 +u - 3 

U.O_ 2 


0.4 ±0.1 


0.69 


0.65 


o.s±8i 


0.2 ±0.1 


0.69 


0.67 


1.4183 


1.4 ±0.1 


0.69 


0.69 


o.7l8:i 


0.6 ±0.1 


0.71 


0.65 


1 n +0 - 5 


1.0 ±0.1 


0.71 


0.67 


i q+0.6 
L -°-0.2 


1.2 ±0.1 


0.71 


0.69 


1-6183 


1.2 ±0.1 


0.71 


0.71 


o.9l8;| 


0.7 ±0.1 


0.73 


0.63 


2 +ai 


0.3 ±0.1 


0.73 


0.65 


i q+0.6 


1.3 ±0.1 


0.73 


0.67 


1 4+0-6 


1.2 ±0.1 


0.73 


0.69 


1 4+0-8 


1.4 ±0.1 


0.73 


0.71 


i -°-0.3 


1.6 ±0.1 


0.73 


0.73 


o.7±8:l 


0.5 ±0.1 


0.75 


0.63 


o.9l8:l 


0.9 ±0.1 


0.75 


0.65 


1-318:5 


1.3 ±0.1 



Table 1: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.75 


0.67 


i o+U.7 
i - z -0.3 


1.4 ±0.1 


0.75 


0.69 


i 7+O.8 
l - 1 -0.2 


1.3 ±0.1 


0.75 


0.71 


o 1 +0.8 
Z - L -0A 


1.4 ±0.1 


0.75 


0.73 


2.018:1 


1.6 ±0.1 


0.75 


0.75 


o.s±8:f 


0.9 ±0.1 


0.77 


0.61 




0.3 ±0.1 


0.77 


0.63 




1.4 ±0.1 


0.77 


0.65 


i.7±8:l 


1.5 ±0.1 


0.77 


0.67 


i c;+0.6 

i -°-0.2 


1.4 ±0.1 


0.77 


0.69 


i 7+O.8 
' -0.4 


1.3 ±0.1 


0.77 


0.71 


l-5la2 


1.4 ±0.1 


0.77 


0.73 


i c+0.8 
i -°-0.3 


1.6 ±0.1 


0.77 


0.75 


2.018:1 


1.6 ±0.1 


0.77 


0.77 


i-i±8:l 


0.9 ±0.1 


0.79 


0.61 


i - u -0.2 


1.0 ±0.1 


0.79 


0.63 


1 fi +0.8 


1.5 ±0.1 



Table 2: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.79 


0.65 


1 7+U.8 
l - 1 -0.4 


1.3 ±0.1 


0.79 


0.67 


i c;+0.8 
i - d -0.3 


1.6 ±0.1 


0.79 


0.69 


1 8+ - 7 
L -°-0.2 


1.4 ±0.1 


0.79 


0.71 


1 8+ - 8 
t-O-OS 


1.5 ±0.1 


0.79 


0.73 


i 7+O.6 
x - ' -0.3 


1.7 ±0.2 


0.79 


0.75 


i -°-0.5 


1.8 ±0.2 


0.79 


0.77 


1 8+ 1 - 2 
i -°-0.5 


1.7 ±0.1 


0.79 


0.79 


i.o±8:l 


1.0 ±0.1 


0.81 


0.59 


0.41°;? 


0.3 ±0.1 


0.81 


0.61 


1-418:5 


1.6 ±0.1 


0.81 


0.63 




1.6 ±0.1 


0.81 


0.65 


1-818:1 


1.5 ±0.1 


0.81 


0.67 


2.018:1 


1.9 ±0.2 


0.81 


0.69 


1 c;+0.8 


1.7 ±0.1 


0.81 


0.71 


i 7+0.9 
x - ' -0.4 


1.5 ±0.1 


0.81 


0.73 




1.9 ±0.2 



Table 3: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.81 


0.75 


9 9 +U.8 
z - z -0.4 


1.6 ±0.1 


0.81 


0.77 


z - 1 - 0.3 


2.2 ±0.2 


0.81 


0.79 


o c+1-3 
z -°-0.5 


2.1 ±0.2 


0.81 


0.81 


1 1+O.6 
L - L -0.2 


1.1 ±0.1 


0.83 


0.59 


i q+0.9 
i -' 3 -0.2 


1.1 ±0.1 


0.83 


0.61 


i -°-0.4 


1.7 ±0.1 


0.83 


0.63 


i 0+O.6 
L -°-0.3 


1.7 ±0.2 


0.83 


0.65 


i c;+0.5 
i -°-0.2 


1.6 ±0.1 


0.83 


0.67 


1 6 +1 '° 


1.7 ±0.2 


0.83 


0.69 


i O+0.9 
J--O-0.4 


1.9 ±0.2 


0.83 


0.71 


z - 1 -0.4 


1.7 ±0.2 


0.83 


0.73 


o i +1.0 
z - 1 - 0.3 


2.0 ±0.2 


0.83 


0.75 


9 9+O.9 
z - z -0.3 


2.0 ±0.2 


0.83 


0.77 


2 3 +L1 

z -°-0.2 


2.0 ±0.2 


0.83 


0.79 


9 9+O.8 
z - z -0.4 


2.3 ±0.2 


0.83 


0.81 


2 6+ 1 ' 2 
z -°-0.3 


2.4 ±0.2 



Table 4: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.83 


0.83 


i q+0.6 
i -°-0.2 


1.3 ±0.1 


0.85 


0.57 


0.3±°;f 


0.3 ±0.1 


0.85 


0.59 


o 9 +i.i 


2.1 ±0.2 


0.85 


0.61 


l-°-0.5 


1.5 ±0.1 


0.85 


0.63 




1.9 ±0.2 


0.85 


0.65 


l-°-0.3 


1.7 ±0.1 


0.85 


0.67 


o o+l.O 
z - z -0.3 


1.8 ±0.2 


0.85 


0.69 


9 n+0.9 
z - u -0.3 


1.7 ±0.2 


0.85 


0.71 


2 4 +o.9 

z - 4 -0.3 


1.9 ±0.2 


0.85 


0.73 


9 q+1.0 


2.1 ±0.2 


0.85 


0.75 


0.5 


2.3 ±0.2 


0.85 


0.77 


9 9+1.1 
z - z -0.5 


2.3 ±0.2 


0.85 


0.79 


2-3lo!3 


2.4 ±0.2 


0.85 


0.81 


z -°-0.4 


2.8 ±0.2 


0.85 


0.83 


o n+1.4 
<J- u -0.5 


3.0 ±0.2 


0.85 


0.85 


1 1+0.8 
1-1-0.2 


1.5 ±0.1 



Table 5: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.87 


0.57 


i q+U.8 
i -°-0.2 


1.4 ±0.1 


0.87 


0.59 


9 n +L1 

Z - U -0.4 


1.8 ±0.2 


0.87 


0.61 




1.7 ±0.2 


0.87 


0.63 


9 n+ L4 


1.8 ±0.2 


0.87 


0.65 




1.9 ±0.2 


0.87 


0.67 


i - y -0.6 


2.0 ±0.2 


0.87 


0.69 


o 9+1-0 
z - z -0.4 


2.2 ±0.2 


0.87 


0.71 


9 9+1.3 
z - z -0.5 


2.1 ±0.2 


0.87 


0.73 


9 fi +0 ' 9 
z - D -0.4 


2.2 ±0.2 


0.87 


0.75 


9 fi +1 - 3 


2.6 ±0.2 


0.87 


0.77 


z -°-0.6 


2.5 ±0.2 


0.87 


0.79 


z -°-0.4 


2.6 ±0.2 


0.87 


0.81 


2.9± ;I 


3.5 ±0.2 


0.87 


0.83 


3.01^ 


3.5 ±0.2 


0.87 


0.85 


q 9+1.4 
°- z -0A 


3.1 ±0.2 


0.87 


0.87 


1 c;+0.6 
i -°-0.2 


1.9 ±0.2 



Table 6: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 



18 



x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.89 


0.55 


0.61^ 


0.4 ±0.1 


0.89 


0.57 


9 9+O.9 
z - z -0.3 


1.9 ±0.2 


0.89 


0.59 




2.1 ±0.2 


0.89 


0.61 


9 q+1.3 


1.8 ±0.2 


0.89 


0.63 


9 c+1.3 
z -°-0.4 


2.1 ±0.2 


0.89 


0.65 


o 9+1.1 

z - z -0.6 


2.3 ±0.2 


0.89 


0.67 


9 9+1.4 
z - z -0.4 


2.3 ±0.2 


0.89 


0.69 


2.6l :i 


2.4 ±0.2 


0.89 


0.71 


9 0+1.2 
^•0-0.4 


2.2 ±0.2 


0.89 


0.73 


9 r+1.4 
z -°-0.6 


2.3 ±0.2 


0.89 


0.75 


O 1+1.1 

« 5 - ± - 0.4 


2.6 ±0.2 


0.89 


0.77 


q 9+1.3 


2.9 ±0.2 


0.89 


0.79 


q 9+1.2 


3.0 ±0.2 


0.89 


0.81 


3 +L2 

O.U_ 4 


3.3 ±0.2 


0.89 


0.83 


q q+1-6 
°-°-0.6 


3.6 ±0.2 


0.89 


0.85 


o fi +1.2 


3.8 ±0.2 



Table 7: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.89 


0.87 


<J- 4 -0.5 


4.7 ±0.2 


0.89 


0.89 


1 q+ - 8 


2.5 ±0.2 


0.91 


0.55 




1.4 ±0.1 


0.91 


0.57 


9 n +1 - 2 


2.0 ±0.2 


0.91 


0.59 


z - 4 -0.3 


2.3 ±0.2 


0.91 


0.61 


9 fi +L4 
z -°-0.5 


2.2 ±0.2 


0.91 


0.63 


2-4± ;I 


2.1 ±0.2 


0.91 


0.65 




2.2 ±0.2 


0.91 


0.67 


9 fi +1 - 5 

Z - D -0.4 


2.5 ±0.2 


0.91 


0.69 


9 ^ +1 - 3 

z -°-0.5 


2.7 ±0.2 


0.91 


0.71 


o 7+1.3 
z - ' -0.6 


2.3 ±0.2 


0.91 


0.73 


3.ll i 


3.0 ±0.2 


0.91 


0.75 


3.0± ;I 


3.3 ±0.2 


0.91 


0.77 


3 c i+ 1 - 8 


3.2 ±0.2 


0.91 


0.79 


3 6 +1 ' 9 


3.6 ±0.2 


0.91 


0.81 


q q+1.6 
O.O_ .4 


3.3 ±0.2 



Table 8: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


x 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.91 


0.83 


q q+1.2 


3.9 ±0.2 


0.91 


0.85 


O + 1.5 

-J-O-0.5 


5.0 ±0.3 


0.91 


0.87 


4 6+ 1 - 4 


5.0 ±0.3 


0.91 


0.89 


4 7+ 1 - 4 

4 - ' -0.5 


5.4 ±0.3 


0.91 


0.91 


o 9+0.7 
z - z -0.3 


2.8 ±0.2 


0.93 


0.53 






0.93 


0.55 


2.7; ; 5 


2.7 ±0.2 


0.93 


0.57 




2.2 ±0.2 


0.93 


0.59 


9 0+1.3 


2.2 ±0.2 


0.93 


0.61 


9 7+I.8 
z - ' -0.7 


2.3 ±0.2 


0.93 


0.63 




2.4 ±0.2 


0.93 


0.65 


2.e; i 


2.3 ±0.2 


0.93 


0.67 


3-0lo1 


2.8 ±0.2 


0.93 


0.69 


3.ol ;i 


3.0 ±0.2 


0.93 


0.71 


3 4 +L7 


3.2 ±0.2 


0.93 


0.73 


3 4 +L1 

J - 4 -0.3 


3.5 ±0.2 



Table 9: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . The 
prediction for the X 3 = 0.93, X 4 = 0.53 bin, for which the theoretical statistical 
uncertainty exceeds 20%, is not included. 
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x 3 


x 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.93 


0.75 


q q+1.6 
°-°-0.6 


3.3 ±0.2 


0.93 


0.77 


O + 1.8 

O.0_o.4 


3.7 ±0.2 


0.93 


0.79 


4 1+ 1 - 7 

0.5 


3.6 ±0.2 


0.93 


0.81 


3.4; ; 5 


4.2 ±0.2 


0.93 


0.83 


3.9lfcg 


5.1 ±0.3 


0.93 


0.85 


4-4± ;l 


5.5 ±0.3 


0.93 


0.87 


4.4± J 


5.6 ±0.3 


0.93 


0.89 


5-ol :i 


6.4 ±0.3 


0.93 


0.91 


4.9+ ; 3 7 


7.7 ±0.3 


0.93 


0.93 


o O+0.9 


3.2 ±0.2 


0.95 


0.53 


1 9+O.6 




0.95 


0.55 


i q+0-9 


1.5 ±0.1 


0.95 


0.57 


2-Hoi 


1.5 ±0.1 


0.95 


0.59 


2.0± ;I 


1.4 ±0.1 


0.95 


0.61 




1.7 ±0.1 


0.95 


0.63 


i - y -0.5 


1.9 ±0.2 



Table 10: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . The 
prediction for the X 3 = 0.95, X 4 = 0.53 bin, for which the theoretical statistical 
uncertainty exceeds 20%, is not included. 
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x 3 


X 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.95 


0.65 


o q+1.3 
z -°-0.3 


1.8 ±0.2 


0.95 


0.67 


9 9+1-2 


2.3 ±0.2 


0.95 


0.69 


o c+1.0 
z - d -0.3 


2.6 ±0.2 


0.95 


0.71 


o 1+1.6 
^■ l -0A 


2.5 ±0.2 


0.95 


0.73 


o 0+1.4 
°-°-0.5 


3.2 ±0.2 


0.95 


0.75 


o 7 +1.9 
°- ' -0.5 


3.9 ±0.2 


0.95 


0.77 


o 0+I.8 
d -°-0.6 


3.8 ±0.2 


0.95 


0.79 


4.ol :l 


5.2 ±0.3 


0.95 


0.81 


4.9;^ 


5.4 ±0.3 


0.95 


0.83 


5.4±S:g 


5.9 ±0.3 


0.95 


0.85 




6.2 ±0.3 


0.95 


0.87 


r 7+I.6 
°- ' -0.8 


7.5 ±0.3 


0.95 


0.89 


r. c+1.8 
°-°-0.9 


7.9 ±0.3 


0.95 


0.91 


r. 4+1.4 


8.2 ±0.3 


0.95 


0.93 


r 1+1.6 
d - i -0.8 


7.1 ±0.3 


0.95 


0.95 


i q+0.8 


2.0 ±0.2 



Table 11: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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x 3 


x 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.97 


0.51 


o.2±8:i 


0.4 ±0.1 


0.97 


0.53 


0.7±o1 




0.97 


0.55 


0.9181 


0.3 ±0.1 


0.97 


0.57 


0.818:1 


0.4 ±0.1 


0.97 


0.59 


i.il83 


0.3 ±0.1 


0.97 


0.61 


o.9±8:I 


0.3 ±0.1 


0.97 


0.63 


i.oi8:S 


0.9 ±0.1 


0.97 


0.65 


i 1+0.8 
i - i -0.3 


0.6 ±0.1 


0.97 


0.67 


i n +0 - 5 


0.8 ±0.1 


0.97 


0.69 


1 4+0-7 


0.7 ±0.1 


0.97 


0.71 


i c+0.9 
i -°-0.4 


0.8 ±0.1 


0.97 


0.73 


1-518:1 


0.9 ±0.1 


0.97 


0.75 


1-91S3 


1.5 ±0.1 


0.97 


0.77 


1.9lfc§ 


1.4 ±0.1 


0.97 


0.79 


q+1.0 
z -°-0.3 


2.2 ±0.2 


0.97 


0.81 


9 S+ 1 - 2 
2 -°-0.3 


3.4 ±0.2 



Table 12: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . The 
prediction for the X 3 = 0.97, X 4 = 0.53 bin, for which the theoretical statistical 
uncertainty exceeds 20%, is not included. 
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x 3 


x 4 


Measured Cross 
Section (pb) 


NLO Cross Section 
CTEQ4M (pb) 


0.97 


0.83 


q o+l.O 


4.0 ±0.2 


0.97 


0.85 


q q+i-3 


4.6 ±0.2 


0.97 


0.87 


4.33:2 


6.5 ±0.3 


0.97 


0.89 


5.4_J: 2 7 


7.4 ±0.3 


0.97 


0.91 


4 q+i-5 

4 - y -0.9 


7.7 ±0.3 


0.97 


0.93 


q + 1.2 

-J-O-0.9 


5.7 ±0.3 


0.97 


0.95 


9 1+1.1 
zi - 0.7 


3.0 ±0.2 


0.97 


0.97 


0.318:? 




0.99 


0.51 


o.ol8:o 




0.99 


0.53 


o.ol8:o 




0.99 


0.55 


0.118:1 




0.99 


0.57 


1 +0 '° 




0.99 


0.59 


n i + - 1 

u-J—o.o 




0.99 


0.61 


n 1 + - 1 

u-J—o.o 




0.99 


0.63 


U- 1 - 0.1 




0.99 


0.65 


U.-L-O.l 





Table 13: The measured and predicted three-jet production cross section in every 
kinematically allowed bin in the Dalitz plane as a function of X 3 and X 4 . Predictions 
for bins with X 3 > 0.98 are not included. 



25 





X 4 


Measured Cross 
Section (pb) 


0.99 


0.67 


u - 1 -0.1 


0.99 


0.69 


o.il8:l 


0.99 


0.71 


0.318:? 


0.99 


0.73 


0.2I8:? 


0.99 


0.75 


o.3l8:l 


0.99 


0.77 


0.318:? 


0.99 


0.79 


0.2I8:? 


0.99 


0.81 


o.6±8:l 


0.99 


0.83 


o.5l8:? 


0.99 


0.85 


o.5l8:^ 


0.99 


0.87 


0-6181 


0.99 


0.89 


0.918:? 


0.99 


0.91 


1 i+0.3 


0.99 


0.93 


1 1 +0 - 4 
L - L -0.3 


0.99 


0.95 


U.O_ 2 


0.99 


0.97 


u - i -0.1 



Table 14: The measured three-jet production cross section in every kinematically 
allowed bin in the Dalitz plane as a function of X 3 and X 4 . 
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In conclusion, we have presented the first comparison of three-jet event cross 
section variation across the Dalitz plane with predictions from a complete NLO QCD 
calculation. The total cross section is found to be 466 ± 3(stat.)l7o 7 (syst.). The 
data agree in absolute magnitude with theory and with our previous measurements. 
The shape of the theoretical and experimental distributions in the Dalitz plane differ 
somewhat; the difference may give an indication of the size of higher order corrections. 
It appears, for example, that up to NLO the theory predicts more soft radiation than 
the data have in the region where the primary partons are approximately back-to- 
back. The data, especially in the region above X3 = 0.98 where a perturbative 
expansion is not reliable, may be useful input to theoretical models of gluon-emission 
processes. 
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